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Abstract 

The higher dimensional spherical symmetric scalar field collapse problem is stud- 
ied in the light of the critical behavior in black hole formation. To make the analysis 
tractable, the self similarity is also imposed. By giving a new view to the self-similar 
scalar field collapse problem, we give the general formula for the critical exponents in 
higher dimensions. In the process, the explanation of the universality of the critical 
phenomena is given within the self-similar context. 



1 Introduction 



Recently, Choptuik analyzed the spherically symmetric scalar field collapse problem nu- 
merically and found several intriguing facts [Q] . Once the initial data which is characterized 
by one parameter p is prepared, one can calculate the evolution of that initial data. Then 
the problem is the fate of the spacetime. If the amplitude of the field is not so large, the 
imploding scalar wave will scatter away and the flat spacetime remains. On the contrary, 
if the amplitude of the field is large enough, the imploding scalar wave will collapse and 
the apparent horizon is developed. The final state is nothing but a black hole. Let us 
consider the verge of two CclSGS, clS the initial data is characterized by the one parameter, 
there is a critical value p* which separates two regions, the corresponding solution is now 
called "choptuon". The findings by Choptuik are echoing, scaling, and universality. The 
echoing is the feature of the choptuon; which shows a discrete self-similarity 



where p and r are logarithms of proper radius r and central proper time t. Here A ~ 3.4. 
More interestingly, in the super-critical case p > p*, there is the scaling law in the mass 
formula if the initial data close to the critical value enough, i.e., 



where the critical exponent (3 is about 0.37. The universality states that the exponent (3 
and the echoing parameter A are independent of any choice of parameter p of the initial 
data. 

The above critical phenomena are also interesting in the light of the cosmic censorship 
conjecture. Indeed, the motivation of Choptuik lied on the cosmic censorship conjecture. 
The choptuon is characterized by a horizonless (massless black hole) spacetime and the 
occurrence of a strong curvature singularity from the regular initial data. The choptuon 
offers a new numerical counterexample to the cosmic censorship conjecture. 

It may also be related to the quantum evaporating black hole. As is well known, when 
we consider the quantum effects, the black hole is no longer black. It emits the radiation 
whose strength is proportional to the surface gravity. Then, as the evaporation proceeds, 
the rate of the emission of the energy is accelerated. The final state might be choptuon, 
and the scaling law might have some physical significance. 




(1) 




(2) 
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Soon after his discovery, Abrahams and Evans numerically investigated the axially- 
symmetric collapse of the gravitational wave and found that the exponent becomes (3 = 

0. 37 and the echo is characterized by A = 0.6. Furthermore, Evans and Coleman || 
studied the radiation fluid system numerically and obtained (3 = 0.36. In this case, the 
choptuon shows the continuous self-similarity. Thus, they suggested the perturbative 
analysis around this self-similar solution to obtain the exponent. Koike et. al. [|j] realized 
this program and confirmed the results of Evans and Coleman. From these results, one 
can observe that the echo strongly depends on the matter content and hence feature of the 
matter. However, it is expected that the scaling is independent on the matter content, 

1. e., there may be a universality in the strong sense. So, we are concentrating on the 
scaling behavior in this paper. Of course, there is a work by Maison |J which studied the 
general perfect fluid by using perturbation method and found (3 depends on the equation 
of state. This result should be confirmed by numerical calculation. Chiba and one of the 
author (J.S.) |J have also investigated the dependence of the phenomena on the theory, 
more precisely, they analyzed Brans-Dicke theory and concluded that (3 depends on the 
Brans-Dicke parameter and A is not sensitive on that parameter. Even if the universality 
in the strong sense does not exist, the universality in the primary sense is interesting 
enough to investigate and we believe it reflects the nonlinear dynamics intrinsic to the 
gravity. 

Almost of the works mentioned above are numerical works, however, it is important 
to give analytical models to understand underlying physics. At present, we know Roberts 
solution as such a model ]7| . This model shows the critical phenomena in the sense which 
will be explained later. One of the purpose of this paper is to give a new look at to this 
model. Then we can understand the origin of the universality properly. 

As we expect the superstring theory is the theory of everything, it may be meaningful 
to consider higher dimensional collapse problem. Moreover, from the universality point of 
view, the various critical behavior should be examined. Therefore, we shall consider the 
higher dimensional scalar field collapse problem with restriction to the spherical symmetry 
and the self-similarity. If we can find the exact solution in general, it is easy to calculate the 
exponents. However, it seems difficult to obtain the exact solution in higher dimensions. 
Thus, the aim of this paper is to give an alternative view to critical behavior exhibited by 
Roberts solution and as an application to calculate the critical exponents of n-dimensional 
self-similar scalar field collapse. 
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The organization of this paper is as follows. In sec. 2, the Roberts solution is reviewed 
and the meaning of the critical behavior is specified. In sec. 3, we present the basic 
equations in the appropriate form for our purpose. In sec. 4, the origin of the universality 
is explained and the general formula for the critical exponents is given. The final section 
devotes to the summary of our work. 

2 Roberts solution and its critical behavior 

In this section, we review the Roberts solution to make the meaning of the critical 
phenomena precise. We will consider the spherical symmetric collapse of a scalar field. 
Let us take the line element as 

ds 2 = -e 2p dudv + r 2 dn 2 , (3) 

where p and r depend only on u and v. Here, we take dfl 2 as a line element of 2-dimensional 
unit sphere. Then, from the Lagrangian 

S = J d 4 x^g[R - ~<Tc> M <R0] , (4) 

we obtain Einstein equations: 

dudy + ^e 2 ? = 0, (5) 

dlr - 2d u pd u r = -ir(9 u 0) 2 , (6) 

d 2 r - 2d v pd v r = -- A r{d v <P) 2 , (7) 

d u d v p - -^d u rd v r = ^e 2p - ^d u (f>d v (f) , (8) 

and the scalar field equation: 

2d u d v <p + d u {\ogr 2 )d v <P + d v {\ogr 2 )d u <P = . (9) 

It should be noticed that the generalization to n massless scalar fields does not alter the 
structure of the equations. Hence, the critical exponent which we shall calculate is the 
same with the case of n massless scalar fields. 
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Now let us impose a self-similarity on the fields such that 



(— )' = 

V 




r+^Af-vf] = 




f"-2p'f = 




p' + VP" + y 2 U-vf) = 


1 e 2 ? 


(/V)' = 


, 



1^/2 



r = vf(r)), (fO) 
P = p(v), (11) 

= (12) 

where rj — u/v. The above ansatz give equations: 

(13) 

(14) 

(15) 

(16) 
(17) 

where the prime denotes the derivative with respective to rj. Eqs. (14) and (15) yield 

p' = . (18) 

It should be noticed that the solution of eqs. (13) and (18) are independent of the number 
of the scalar fields. Using eq.(18), eq.(15) becomes 

f" = ~\f^ ■ (19) 

Eqs. (13), (18) and (19) give a closed equation system. Eqs. (16) and (17) are redundant 
equations. Without loosing generality, we can put e 2p = 1. Using eq.(13), we obtain 

f = ±{ r] 2 -2r ] +l-p 2 } , (20) 

where p is the integration constant. Substituting the solution (20) into eq.(19), we get 

= ±log T r - 21 

u — (1 — p)v 

Consequently, the line element becomes 

ds 2 = -dudv + ^(u 2 - 2uv + (1 - p 2 )v 2 )dQ 2 . (22) 



Now, we shall see how the critical phenomena arises. It is straightforward to calculate 
the spacetime curvature as 

„ p 2 uv , . 

* = (23) 
The curvature singularity occurs at u — (1 — p)v (v > 0) and u — (1 + p)v (v < 0). The 
apparent horizon is given by the equation 

^'^^-0, (24) 

where 9 is the expansion of the outgoing null rays. Hence, the apparent horizon is located 
at u — (1 —p 2 )v. The structure of the spacetime depends on the parameter p. In the case, 
< p < 1, the singularity becomes timelike in any region. There is no apparent horizon 
and therefore no black hole. Physically, this case corresponds to the weak implosion of 
the scalar wave. In the case, p > 1, the singularity becomes spacelike in the region v > 0, 
while the singularity in the region v < remains to be timelike. Moreover the spacelike 
singularity is hidden by the apparent horizon. Therefore, there is the critical phenomena 
in this system. The critical parameter of this system is p* — 1. To investigate the scaling 
behavior, it is necessary to calculate the mass of the black hole in the supercritical region. 
Due to the spherical symmetry, there is a local mass function 

M= r -(l-g ah d a rd b r) , (25) 

where a, b runs only v and u. On the apparent horizon, this mass (25) becomes 

M,. H . = ?y^L . (26) 

Unfortunately, this quantity diverges in the asymptotic future due to the self-similarity. 
This fact is also understood from the eq.(21) which shows the imploding scalar wave also 
diverges in the asymptotic future. In other words, 

2 

V V 

M u ^-oo = P — (27) 

also diverges. Therefore, we will define the order parameter as 

M A . H . 



M u ^_oo p 



(28) 
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Notice that, in the present case, p* = 1, hence, it shows the scaling behavior 

(P-Pl 1/2 , (29) 



M A . H . 



Mi— oc 

here the exponent can be read off as 1/2 || ||. 

In this sense, we say the self-similar spherical collapse shows the critical behavior 
whose critical exponent is 1/2. In the following discussion on the higher dimensional 
critical phenomena, this definition of the order parameter is assumed. Of course, it is 
possible to make the spacetime asymptotic flat by cutting and pasting the spacetime 
appropriately as is discussed in j|, [J. As we are interested in the critical behavior, not 
the cosmic censorship conjecture, we are not concerned with it anymore. 

3 Collapse dynamics as an autonomous system 

As it seems difficult to find exact solution in higher dimensions by the method similar 
to that used in the previous section, an alternative method must be found to calculate 
the critical exponent in black hole formation. For that purpose, we will use the similar 
formalism used by Brady ]10[ in the slightly different context. Using advanced Bondi 



coordinates v, r, 9i, • • • n _ 2 , the spherical symmetric metric can be written 

ds 2 = —g(v, r)g(v, r)dv 2 + 2g(v, r)dvdr + r 2 dQ 2 l _ 2 , (30) 

where dfl 2 l _ 2 is the line element of the hyper sphere in n-dimensions. Hence r is the 
proper circumferential radius. Using this parameterization, the equations of motion can 
be derived from the action (4). From G rr = 8irT rr , we obtain 

(log£), r = ^-r(0 ir ) 2 . (31) 



n - 2 



The equation G r v = 8ttT^ leads to 



Furthermore, one can find 



+ 9<!>M ■ (32) 



Roe = 8n[T ee - -±—g ee ( g ^T^ + (n - 2)g ee T ee )} = , (33) 
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where we have used the fact that traceless condition of the energy momentum tensor in 
2-dimensional scalar field. From this result, we obtain 

(rg) )r = {n-3)g-{n-4)g . (34) 

And the equation of motion of the scalar field A0 = gives 

[gr n - 2 (t>A,r = -in - 2)r"- 3 it , - 2r^ 2 0,™ • (35) 

The remaining equations can be derived from the above equations. These equations 
motion are still difficult to solve, so we further restrict the model to be self-similar. More 
precisely, we impose the condition 

L^g^v = 2g iXV , (36) 
where £ = rd r + vd v . Now, define x = r/v , then we obtain the function forms as 

f x i(v) 

g = g{x), g = g(x), <f> = — dn. (37) 

Jo n 

In general, the scalar field has more freedom, however, we do not consider it in this paper. 
Indeed, that case has not do with the critical behavior we have mentioned in sec. 2. Using 
the above results, equations of motion can be written as 

xg' = ^2#7 2 , (38) 

g> = ^Zl^-g) (39) 

x 

9-9 = — f 1? -{g-2x) 1 \ (40) 

[n — 2)(n — 3) 

x(g-2x)i = -j(n-3)g- (n-2)x , (41) 
where a prime denotes the derivative with respect to x. Following Brady HlOl , we define 



V = - , z= - , £= -logx . (42) 
9 9 

Here we have taken into account the fact that when v varies from to oo, x varies from 
oo to and hence £ varies from — oo to oo. Thus, the basic equations are deduced as 

y = _( n _3) + [-^ 7 2 + (n-3)]y , (43) 



z = -z[n - 2 - ^— -] , (44) 
(l-2z)<y = 1 [ V —^-{n-2)z], (45) 



7 = ±< 



(n-2)(n-3)(l- y: ^) ^ 



8tt(1-2*) 

where a dot denotes the derivative with respect to £. Due to the existence of the con- 
straint equation, the apparently 3-dimensional system reduces to 2-dimensional dynamical 
system. 

It is easy to find the fixed point of this system 



n — 3 n — 2 

Vf = ^=2> */ = 1 > ^ = ± V^T- (47) 

In Fig.l, we have shown the phase flow in the 4-dimensional case. It certainly shows the 
fixed point is a saddle point. And the apparent horizon is given by 

g^dprdvT = - = y = . (48) 

Then, at the fixed point, the phase flow bifurcate to the black hole and the flat spacetime. 
The line which separate both type might be called as choptuon. As for the comparison, 
the result of 6-dimensional case is also presented in Fig. 2. It shows almost the same 
behavior with that of 4-dimensions. 

4 Universality of the critical phenomena 

As the results of the previous section indicate, the phase flow structure is determined by 
the fixed point. The important point here is that the fixed point is the saddle point, hence 
the relevant mode explained below is unique. This explains the universality of the critical 
behavior in black hole formation. It should be noted here the similarity of our results 
with those of Koike et. al. |IJ] and Maison |5|]. However, in our case as opposed to theirs, 
the global phase structure is known. So explanation is complete. In a sense, the model 
considered in this paper is the minisuperspace model for critical phenomena. 

To make the concept clear, it is convenient to define several important operator. The 
time-evolution operator is defined as 

RtX(0) = X(0 , (49) 
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where X = (y, z). Then the fixed point is characterized by 



R^Xf - X f . 



(50) 



The tangent map of at a fixed point Xf is defined as 



„ , RclX f + eF) -X f 
F = lim €V J ^ ^ 



(51) 



An eigenmode Fj of can be defined as 



(52) 



A mode with k > is called as a relevant mode from an apparent reason. In the case of 
k = 0, we say it is marginal and if k < we say irrelevant. Now what we have shown 
in the previous section is the existence of the saddle fixed point and the critical line on 
which the every point goes to the fixed point. Hence there is a one relevant mode as is 
mentioned. 

First of all, we must prepare the initial data, /, characterized by the one parameter 
p. A one parameter family of initial data / will in general intersect with the critical line, 
and the intersection point, X c will be driven to the fixed point under the time evolution: 



Now, we can describe why the universality exists in black hole formation (see Fig.3). The 
intersection point is nothing but the critical point, corresponding to the parameter p*. So 
the critical line corresponds to the critical solution. If we start from the point which is 
different from the critical point but close to it, it will approach to the fixed point along 
the critical line. However, at near the fixed point, it will eventually be driven away along 
the direction corresponding to the relevant mode. If the initial parameter is larger than 
the critical value p*, it goes to the black hole and otherwise it goes to the flat spacetime. 
Apparently this behavior is determined only by the behavior around the fixed point, hence 
it is independent on the choice of the initial parameter. 

Using this picture, we can give the general formula for the critical exponents in n- 
dimensional spacetime. Let X init be the initial data such that 




. 



(53) 



Ai n i t — X c — eF , 



(54) 
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where F is the vector which deviate the initial data from the critical value, and e =| p—p* |. 
What we are interested in is the time-evolved value of this initial data, 

i^Xinit = R^Xc-eF) (55) 
= X c (0-eT 5 F + O(e 2 ) (56) 
= Xf-ee^F^, (57) 

where we used eqs. (51) and (53) in the above manipulation. In the final line, the relevant 
mode F rel becomes dominant. Thus we have obtained the asymptotic behavior of the time 
evolved initial data. Let us recall the apparent horizon is given by y = 0, hence 

y f = ee*F iel , y (58) 
= 6(^rF rel , y , (59) 

where we have used the relation £ = — \ogr/v. Here, F re i iy denotes the y-component of 
the relevant mode F rei . Therefore, we obtain the scaling relation 

t a.h. ~ • (60) 
From the dimensional consideration, the black hole mass scales as 

M b .h. ~ rTl ~ ^ ~l V - P* f , (61) 

where we have identified as j3 = (n — 3)//c. The remaining task is the calculation of k 
by using the linear perturbation. By substituting the linear expansion around the fixed 
point 

V = Vf + $y ,z = z f + 5z ,7 2 = 7^ + <57 2 , (62) 

into eqs. (43), (44) and (46), we obtain 

n — 3 

5y = 2^5z, (63) 

5z = { - Uy, (64) 

n — 3 

2 (n-2) 3 2(n-2) 

57 = ~8^(^3) 8^ ' (65) 
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It turns out that the fixed point is saddle point from these equations. It is also straight- 
forward to calculate the eigenvalue of the relevant mode as 



k = yJ2(n - 2) . (66) 

Finally, we obtain 

p(n - 2) 

Thus, we have completed our task. 



5 Conclusion 

We have investigated the higher dimensional spherically symmetric self-similar scalar field 
collapse from the point of view of the critical behavior. In the course of the analysis, we 
found the completely satisfactory explanation of the universality of the critical behavior 
in black hole formation within this restricted context. Namely, the universality is nothing 
but the reflection of the phase flow structure of the system. In the system considered in 
this paper, there is a unique fixed point which has only one relevant mode. Then, the 
asymptotic behavior of the system is determined by the behavior around the fixed point. 
Hence, the critical behavior is independent on the choice of the parameter p. This result 
leads us to the formula for the critical exponents in higher dimensions. From our formula, 
it is easy to calculate the critical exponent as in the Table. 1. In the 4-dimensional case, 
it completely agrees with the result obtained by using the exact solution. 

It is interesting to consider the relation between our work and that of Tomimatsu fTI 



concerning with the quantization of the self-similar spacetime. In the future, we would 
like to attack to this problem. 
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dimension n 


exponent (3 


4 


0.5 


5 


0.81650 


6 


1.06066 


7 


1.26491 


8 


1.44338 


9 


1.60357 


10 


1.75 


11 


1.88561 


12 


2.01246 



Table 1: This table shows critical exponents for several dimensions. In the 4-dimensional 
case, the exponent (0.5) completely agrees with the result obtained by using the exact 
solution. In higher dimension, the exponent takes larger value. Especially n > 6, the 
exponent is larger than 1, this fact shows that black hole mass is smoothly connected to 
zero when p — > p*. 
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Figure captions 



Fig.l: Phase flow in the 4- dimensional case is presented. Fixed point of the 4- 
dimensional case is at y = 0.5, z — 1. On z — 2, initial data are set. And there, ys 
play the role of initial parameter. The flow curves which start from y smaller than critical 
y bifurcate to the black hole, the others bifurcate to the flat spacetime. 

Fig. 2: Phase flow in the 6-dimensional case is presented. Fixed point of this case is 
y = 0.75, z = 1. Other features are not different from the 4-dimensional case. 

Fig. 3: Schematic figure of phase flow is shown. Solid curves shows the phase flow, 
and start from I which is a one parameter family of initial data. They go along critical 
line and bifurcate at fixed point Xf. 
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Fig. 1 




i 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

y 



N 



2 r 

1.8 - 

1.6 - 

1.4 - 

1.2 - 

1 - 

0.8 - 

0.6 - 

0.4 - 

0.2 - 

- 








0.2 



0.4 



0.6 



0.8 




Fig.3 



